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Abstract 

This paper proposes a parametric approach for stochastic modeling 
of limit order markets. The models are obtained by augmenting clas- 
sical perfectly liquid market models by few additional risk factors that 
describe liquidity properties of the order book. The resulting models are 
easy to calibrate and to analyze using standard techniques for multivari- 
ate stochastic processes. Despite their simplicity, the models are able to 
capture several properties that have been found in microstructural anal- 
ysis of limit order markets. Calibration of a continuous-time three-factor 
model to Copenhagen Stock Exchange data exhibits e.g. mean reversion 
in liquidity as well as the so called crowding out effect which influences 
subsequent mid-price moves. Our dynamic models are well suited also for 
analyzing market resiliency after liquidity shocks. 



1 Introduction 

Most modern stock exchanges are based on the continuous double auction mech- 
anism where outstanding limit orders are organized in the limit order book. If 
the size of a market order exceeds the quantity available at the best price, the 
trader has to climb up the book and accept worse marginal prices to get his 
order filled. In other words, marginal prices of market orders are increasing 
functions of the order size. This is in contrast with classical market models 
where marginal prices are stochastic processes independent of order size. 

This paper presents simple extensions to classical perfectly liquid market 
models in order to describe liquidity aspects of limit order markets. The models 
are obtained by introducing few additional risk factors that describe the liquidity 
properties of the order book. This allows for compact descriptions of both 
immediate liquidity costs of market orders as well as resiliency properties of the 
order book. Our models can incorporate a wide variety of existing models for 
market prices while leaving room for modeling liquidity aspects of the market. 
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Classical perfectly liquid market models are obtained as special cases when 
illiquidity effects vanish. 

There exists an extensive literature on microstructural models where ar- 
riving orders are related to various features of the limit order book; see e.g. 
Glosten [IT], Biais et al. [SJ, Parlour [37], Sandas 132;;, Luckock [35], Smith et 
al. [33] j Ranaldo [3D], Foucalt et al. [TBJ, Biais, Glosten and Spatt [1], Bouchaud 
et al. [6J, Rosu [31] and their references. Such detailed studies provide valuable 
insights to the formation of the order book but the employed models may be 
cumbersome in describing the statistical behavior of the book over time. Re- 
cently, Cont et al. introduced a more tractable model where orders arrive 
at discrete price levels according to independent Poisson processes with a rate 
depending on the distance to the best quotes. Their model gives a stochastic 
description of the complete order book which is essential e.g. when studying 
the costs of trading strategies. Several empirically observed phenomena were 
reproduced by the model of [TT] in numerical calculations employing Laplace 
transform techniques. 

The present paper proposes a further simplification by describing order book 
dynamics through a finite number of risk factors without explicit account of the 
order flow. One of the factors is the market price while the others describe liquid- 
ity properties of the order book. This reduction provides a natural link between 
classical perfectly liquid market models and microstructural models of limit or- 
der markets. With appropriate parameterization, essential features of the order 
book can be captured already with low-dimensional models. The dimensionality 
becomes an important issue, for example, in numerical optimization of trading 
strategies; see e.g. Bcrtsimas and Lo [3] or Koivu and Pennanen [35] . 

Another advantage of our approach is its simplicity. The resulting models 
are easy to calibrate and to analyze using standard techniques for multivariate 
stochastic processes. This is illustrated by calibrating continuous-time order 
book models to limit order data from Copenhagen stock exchange. We model 
the market price and two liquidity factors by a three-dimensional stochastic 
differential equation which is calibrated to historical observations of order books. 
Despite its simple structure, the model is able to reproduce some well-known 
properties of limit order markets. In particular, the calibrated models are mean 
reverting in liquidity and they exhibit the so called "crowding out" effect, which 
means that high liquidity on one side of the book encourages more aggressive 
limit orders and quote improvements on that side. In our reduced form model, 
this shows up as increases in the mid-price when the ask-sidc liquidity, in terms 
of the liquidity factors, is high compared to the liquidity on the bid-side and vice 
versa. Both effects are supported by extensive theoretical and empirical studies 
of limit order markets; see e.g. Biais et al. [SJ, Foucault, Kadan and Kandel [TP] . 
Parlour [37], Griffiths et al. [TS], Ranaldo [3D], Bouchaud et al. [DJ. 

Our models are close to those developed in Cetin, Jarrow and Protter [DJ and 
Cetin, Jarrow, Protter and Warachka [8] . The main difference is that our models 
retain the monotonicity of marginal prices of market orders which implies that 
their total costs are convex. Convexity is an essential feature of limit order mar- 
kets and it has many important implications for risk management and pricing 
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and hedging of financial instruments; see e.g. Delbaen and Schachermayer [13] 
or Follmer and Schied [15] for classical perfectly liquid models and Jouini and 
Kallal [22] or Kabanov [23J for models with proportional transactions costs. Our 
models can be seen as generalizations of the model of Cetin and Rogers [TU] who 
studied utility maximization under liquidity costs given in terms of a fixed con- 
vex function of the order size. Basic arbitrage and pricing theory for general 
convex cost functions (that cover, in particular, piecewise linear functions as- 
sociated with limit order markets) has been developed in Pennanen [251 [55] . 
Besides theoretical aspects, convexity is essential also in optimization of trading 
strategies; see e.g. Edirisinghe, Naik and Uppal [T3], Bertsimas and Lo [3J, Alm- 
gren and Chriss [2 , Huberman and Stanzl [2T] , Alfonsi, Schied and Fruth pQ or 
Koivu and Pennanen [25] . 

The rest of this paper is organized as follows. Section [2] gives a short review 
of the order book structure and a quantification of liquidity costs in monetary 
terms. This is used in Section [3] to describe liquidity costs by simple paramet- 
ric models driven by few liquidity factors. The stochastic order book models 
are then obtained by specifying stochastic models for the behavior of the fac- 
tors. Section |4] presents an empirical study using data from Copenhagen stock 
exchange and Section [5] concludes. 

2 Liquidity costs of market orders 

We begin with a quick review of order book structure in order to introduce 
the main concepts, the data and the notation. The limit order book (LOB) 
maintains a record of all submitted limit orders that have not been cancelled 
or met by a market order. Table [T] displays a part of a limit order book of 
one particular stock at the Copenhagen stock exchange at a given time instant. 
When submitting a buy market order, only a finite number of shares can be 
bought at the best ask price and when buying more, one gets the second lowest 
price and so on. For example, exactly 20800 shares could have been bought at 
the best ask price of 239 in the LOB of Table [1] 

The marginal price of a buy market order is a positive, nondecreasing, piece- 
wise constant function of the number of shares bought. For sell market orders, 
the situation is similar and the marginal revenue of selling is a positive, nonin- 
creasing, piecewise constant function of the number of shares sold. Interpreting 
sell orders as buy orders of negative quantity, we can, as e.g. in Glosten [17] , 
incorporate the bid and ask sides of the book into a single curve s : R M> [0, +oo] 
that gives the marginal price of a market order of arbitrary sign and quantity. 
Since the ask price is always strictly greater than the bid price, s is monotone 
and positive. Figure Q] plots the marginal price curve associated with the LOB 
in Tabled] 

The monotonicity of the marginal price curve s is essential when studying 
trading costs. The positivity and monotonicity of s imply that the total cost 
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Table 1: Example limit order book. The table presents a snapshot of TDC A/S 
limit order book observed at 12JAN2005:13:58:19.430. The book is obtained 
from Copenhagen Stock Exchange order flow data by implementing the rules of 
SAXESS trading protocol. 
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Figure 1: Marginal price curve corresponding to Table [T] The horizontal axis 
gives the cumulative depth of the book measured in the number of shares. 
Negative order quantity corresponds to a sale. The price unit is 1 DKK. 
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of a market order of x shares is a nondecreasing convex function that vanishes 
at x = 0. A negative x incurs a negative cost which just means that sales yield 
revenue. Since s is piecewise constant, S is piecewise linear. In a perfectly 
liquid market, the marginal price s would be constant and the cost function 
S would be linear. Consequences of nonlinearities on arbitrage theory and 
on pricing and hedging have been studied in a number of papers. Jouini and 
Kallal [22 studied the effect of a bid-ask spread with infinite supply at the best 
quotes. Cetin, Jarrow and Protter [5] studied markets with nonlinear illiquidity 
effects in terms of the "supply curve" , which in our notation corresponds to the 
function x >—¥ S{x)/x. Cetin and Rogers [TU] studied portfolio optimization in 
markets with diffcrcntiable convex cost functions. General convex cost functions 
that cover, in particular, piecewise linear functions associated with limit order 
markets have been studied in Pennanen [28j [29] . There the aim was to study 
the effects of nonlinearities on arbitrage and superhedging theories. 

Perfectly liquid market models are often taken as approximations of real 
markets where the marginal price curve is replaced by a horizontal line passing 
through the mid-price ("market price") 

s b + s a 

where s b and s a are the best bid and best ask price, respectively. In terms of 
the cost of market orders, 

s b .- lim and s a . = Um M. 

t/O X x\0 X 

Modeling of s is a classical subject in financial econometrics. In this paper, 
we focus on modeling the nonconstancy of the marginal price curve around the 
mid-price s. 

From practical point of view, it is more interesting to describe liquidity 
in terms of monetary units rather than in numbers of shares. For example, 
one would not expect the shape of an order book to remain the same after a 
stock split; see also [8] Section 4.4]. While market values of shares are usually 
regarded as upward drifting processes it might be that illiquidity effects are of 
more stationary nature in the long run when measured in monetary units. A 
simple monetary measure of illiquidity is given by the function 

r(h) :=ln^-^ =lns(h/s) -Ins 
s 

which gives the percentual change in the marginal price relative to the mid-price, 
as a function of the mark-to-market value h = sx of a market order of x shares. 
Taking logarithms of the whole marginal price curve is analogous to the common 
practice of modeling market price dynamics by describing the behavior of its 
logarithm. 

We will call r the relative price impact curve. It describes the temporary 
price impact that a market order would have on the best quotes (bid or the 
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Figure 2: The price impact curve r (step function) corresponding to Table [T] 
together with its linear approximation. The horizontal axis gives the mark-to- 
market value h of the order quantity. The unit of h is DKK 1 Million. 



ask, depending on the sign of the order). It is always nondecreasing and passes 
through the origin. In perfectly liquid markets r = 0, whereas steeper r corre- 
sponds to illiquidity in execution of market orders. Several authors have studied 
how trades affect the expected mid-price at some point in the future; see e.g. 
Hasbrouck 20 j, Bouchaud et al. [B] Section 5] and the references therein. Such 
dependencies are often called "price impacts" but they should not be confused 
with the relative price impact curve r which describes the directly observable 
instantaneous impacts of market orders. The relative price impact curve r corre- 
sponding to the LOB given in Table Q] is drawn as the step function in Figure [2j 

The marginal price curve s and the cost function S can be expressed in terms 
of the mid-price s and the relative price impact curve r as 



respectively. Making the change of variables x — h/s, we get the representation 



= se 



op r ( sx ) 



and 




S(x) — ip{sx) 



(1) 



where 




The function ip is convex, nondecreasing and vanishes at the origin. It satisfies 
(f(h) > h, where equality holds in perfectly liquid markets. Indeed, since r is a 
nondecreasing function passing through the origin, we have e r ^ h ' > 1 for h > 
and e r< - h ^ < 1 for h < 0. Integrating from to h gives <p{h) > h. Equality holds 
exactly when r = 0. 

The function tp gives the cost of a market order whose mark-to-market value 
is h. The decomposition of S in (QJ supports the common practice of regarding 
market risk and liquidity risk as separate sources of financial risks. Whereas 
unpredictable variations in the mid-price s are often interpreted as market risk, 
variations in <p (or, equivalently, in r) may be interpreted as liquidity risk. In 
the remainder of this paper, we focus on modeling the liquidity cost function ip 
and its relation to the mid-price s. 



3 Reduced form LOB models 

In classical perfectly liquid market models, the mid-price s is a positive stochas- 
tic process and it is assumed that marginal prices of market orders are inde- 
pendent of order size. This corresponds to modeling the limit order book by a 
horizontal line passing through the mid-price s so that r = 0. In this section we 
propose simple extensions where the shape of the limit order book is described 
by modeling the relative price impact curve r as a simple parametric function 
of the order size. The parameters may be interpreted as risk factors that de- 
scribe liquidity properties of the market. A stochastic order book model is then 
obtained by describing the mid-price and the liquidity factors as a multivariate 
stochastic process. 

A minimal condition for a sensible market model is that the marginal price 
curve s is positive and nondecreasing, or equivalently, that the cost function S is 
increasing and convex. This can be achieved by expressing the cost function in 
terms of the relative price impact curve as in (JT|) and modeling the price impact 
curve r as a nondecreasing curve passing through the origin^. A simple model 
that looks beyond the bid and ask prices is obtained with 

r(h) = (3h, (2) 

where (3 > describes the overall degree of illiquidity in the book. Larger the 
value of /3, less liquid the market. This is reminiscent of the linear specification of 
the "price impact function" in Sandas [32]. However, the price impact function 
in Sandas [32] describes the dependence of the "fundamental value" of the stock 

1 This is parallel to the common practice of modeling the term structure of interest rates by 
describing the behavior of the forward curve. The positivity of the forward curve guarantees 
that the zero curve is decreasing as a function of maturity and that it lies between zero and 
one. 
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(one period ahead) on the order quantity while the relative price impact curve 
r describes the shape of the limit order book which is directly observable. 
The solid line in Figure [2] corresponds to the average slope 

f h+ 

P := argmin / \r(h) — (3h\ 2 dh 

/36R Jh~ 

of the relative price impact curve of the LOB in Table [TJ Here h + > and 
h~~ < are finite cut-off parameters that limit the focus on the best 10 prices 
on each side of the book. The linear model ^ corresponds to 

<p(h) = / e^dz = — — 
Jo P 

and thus, to the cost function 

fisx _ 1 

S(X) = 



This is close to the exponential model studied in Cetin and Rogers [TU] ■ In their 
model S(x) = s(e ax — l)/a which describes liquidity effects in terms of units 
of shares and implies, in particular, that iUiquidity effects diminish when the 
market price s increases and vice versa. 

In some situations, the bid and ask sides of a limit order book exhibit dif- 
ferent degrees of illiquidity and the numbers 

f h+ 

f3 + := argmin / \r(h) — fth\ 2 dh, 
pm Jo 

,o 

/3 := argmin / \r{h) - (3h\ 2 dh, 
pm Jh- 

may be significantly different. Differences in liquidity on the two sides of the 
book is known to be significant in explaining e.g. mid-price movements; see 
Parlour [27], Griffiths et al. [18], Ranaldo [30] and Cont et al. [H]- We will 
model the relative price impact curve by the piecewise linear function 

r(h) = \^ +k fOT ^°' (3) 
V ; \p~h for h < 0. W 

Figure |3] plots the corresponding approximation of the LOB in Table [1] This 
corresponds to the cost function 



S(x) 



p + for x > 0, 
for x < 







This could be extended in an obvious way to cover nonzero bid-ask spreads or 
proportional transaction costs as e.g. in Jouini and Kallal [22] (who assumed 
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Figure 3: The price impact curve r and its two sided linear approximation. The 
horizontal axis gives the mark-to-market value h of the order quantity. The unit 
of h is DKK 1 Million. 



infinite supply at the best quotes) . This would correspond to a discontinuity of 
r at the origin. The bid-ask spread could be taken as an additional liquidity 
factor in the model. 

Given a parametric description of the order book, we obtain stochastic mar- 
ket models with illiquidity effects by modeling the parameters as a multivariate 
stochastic process. If the relative price impact curve is parameterized by a 
vector 9 £ K d , then the marginal price curve is given by 

s{x) = se r{Sx ' e) . 

The corresponding cost function can be expressed as in ([1]) as 

S(x) = / e rl - z ^dz. 
Jo 

We can employ the wide variety of existing stochastic models for the mid-price 
s and augment them with stochastic descriptions of 0. The corresponding cost 
functions define a convex cost process in the sense of [25J and are thus amenable 
to analytical techniques similar to those applied to classical perfectly liquid 
market models; see [28] and [29]. The concrete interpretations of the parameters 
are useful when specifying the stochastics. 

Model 1 // the relative price impact curve is parameterized by ([3]) 7 then the 
market model is driven by the three-dimensional process (s, f3~ , f3 + ). A simple 



stochastic specification is to assume a univariate process for the mid-price s and 
that the two dimensional vector 9 = (ln/3 _ , ln/3 + ) follows a two-dimensional 
Ornstein- Uhlenbeck process 

d9 = T( f i-8)dt + Y,dW t , 

where 9 £ l 2x2 , /i £ I 2 and E £R 2x2 are parameters of the model and W is 
e.g. a two-dimensional Brownian motion. 

The following allows for modeling dependencies between the mid-price and 
the liquidity factors. 

Model 2 Assume that the three-dimensional process t; = (In s, ln/3 - , ln/3 + ) sat- 
isfies the linear stochastic differential equation 



d£ = {A£ + a)dt + Y.dW t . 



(4) 



where A £ R Jx3 and a £ K. . In particular, the mid-price drift may depend 
on the shape of the order book through the liquidity factors f3~~ and f3 + ; see 
Section fjj The solution of ^ can be written as 
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+ / e^ t ~ s '> A ads + / e^-^ZdWt 



It seems reasonable to assume that the drift is stationary in the sense that there 
is a S £ R 3 such that 

lim E[A£ t + a] = S. 

t— >oo 

In particular, we would expect that the liquidity factors are stationary so that 
the second and the third component of S are zero. Unless the mid-price is also 
stationary, we must then have that the first column of A is zero since otherwise 
the expectation could not exist. The average values of the liquidity factors are 
then determined by 



A 22 
A 32 



A 23 
A 33 





E/3- 




a 2 









E(3+ 


+ 


as 








and the average mid-price drift becomes 

Si = A 12 Ef3- + A 13 E/3+ 



ai. 



4 Empirical study 

As an illustration of the modeling approach described above, we will calibrate 
the three-dimensional Model [2] to order book data from Copenhagen Stock Ex- 
change. We first adjust the liquidity factors for intra-day patterns that are 
found to explain a considerable part of the variations. The three-dimensional 
stochastic differential equation is then calibrated to discrete observations of the 
mid-price and the deseasonalized liquidity factors. The resulting model exhibits 
certain well-known features of limit order markets. 
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4.1 Data 



We use a data set from Copenhagen Stock Exchange (CSE) that contains the 
buy and sell orders and transactions of all traded stocks for the period January 
2005 to March 20050 Each order record carries information about the date 
and time of submission, the order type, the quantity and the price. Trading 
on CSE takes place in the common trading system SAXESS, which implements 
a continuous double auction mechanism. In this study we focus only on the 
trading lot market and have removed all odd lots from the dataset. For a 
detailed description of the trading conditions and characteristics of the system 
we refer to NOREX Member Rules book[l 

We use the event history to construct the order books and the corresponding 
mid-price s and the bid and ask side liquidity factors /3~ and /3 + as described 
in Section [2] The computations are carried out using SAS and R computing 
environments. Starting from the initial state of the order book at the beginning 
of the trading day, we update the book at every new event by implementing 
the rules of SAXESS trading protocol. All new orders, deletions and updates of 
earlier orders and trade executions are accounted for as they arrive. We use the 
order's direction, price and position to determine its rank in the book according 
to price/time priority. When constructing the books, we do not include hidden 
volume (contained e.g. in iceberg orders) in our computations to ensure that 
the books correspond to the information available to traders over time. Also, 
to eliminate the effects pre- and post-trading sessions, we exclude the first and 
the last trading hours from the data. 

4.2 Intra-day patterns in liquidity 

As discussed e.g. in Engle and Russell [14], the frequency of trades, volume, 
and spreads all typically exhibit consistent patterns over the course of the day. 
Therefore, before calibrating our model, we account for intra-day liquidity pat- 
terns by regressing In (3^~ and In /3 t ~ on trading hour dummies. Since the liquidity 
factors are all strictly positive we have taken logarithms before the linear re- 
gressions. The results for 10 stocks are plotted in Figure [4] All stocks exhibit 
decreasing values of j3 + and /3~ over the trading day. This means that the 
liquidity of the order book tends to increase as more orders are accumulated 
over the day. This is consistent with earlier studies of intra-day patterns in 
order books. For example, Kempf and Mayston [531 Section 2] observed that 
the spread decreases and the depth of the book increases over the trading day. 

Another notable pattern is that the bid side seems steeper on average than 
the ask side. This may reflect the fact that traders on the buy side are exposed 
to market crashes which are often perceived more likely than sudden upward 
moves of similar magnitude. On the other hand, given the relatively short time 

2 We thank Nikolaj Munck at OMX market research in Stockholm for providing the data 
set and assisting with issues involved in reconstruction of the limit order books. 

3 See http: / /www. nasdaqomx.com/listingcenter / nordicmarket / rulesandregulations/copenhagen 
for the latest copy of NOREX Member Rules. The results in this paper are based on the 
rules effective over the data period January-March 2005 
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period, this phenomenon may also be related to market conditions and may be 
reversed e.g. in bull markets. 

From now on, the symbols (3 + and f3~ will denote the deseasonalized vari- 
ables. Figure [5] plots the evolution of the logarithms of /3 + and f3~ together with 
contour plots of the kernel density estimate of their joint probability density for 
TDC's limit order book. The series is obtained with 10 minute sampling inter- 
val. The liquidity factors present significant variation over time. The findings 
are similar for the rest of the stocks. Figure [5] plots the autocovariance functions 
of the liquidity factors for TDC and Moeller-Maersk. The dynamic behavior of 
the factors are characterized by strong positive autocorrelations which decay 
exponentially. 

4.3 Calibration of a linear model 

Consider Model[2]where the order book is driven by the three dimensional vector 
£ = (In s, In /3~, ln/3 + ). For I > 0, the conditional distribution of £z given £o is 
Gaussian with mean 

£ = e<%+ f e^ A ads 
Jo 

and variance (by Ito isometry) 

Vi = f e^-^^-^Yfds. (5) 
Jo 

In other words, 

6 = B^ + bi + e, (6) 

where B\ = e lA , 

bi= [ e^ A ads 
Jo 

and e is Gaussian with zero expectation and variance V/. The model can be 
calibrated to discrete data with fixed time increments I by first estimating the 
values of Bi, bi and Vi and then solving for the corresponding A, a and E. 
Indeed, © can be seen as a linear time series model whose parameters can 
be estimated with ordinary least squares. Once estimates of Bi, bi and Vi are 
available, the matrix A can be estimated with the matrix logarithm 

A = -InA- 

On the other hand 

a = A^bi, 

where the matrix 

Jo 
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Figure 4: Intra-day patterns in the liquidity factors. The two lines in each 
figure give the estimated hourly dummies for the liquidity factors. The solid 
line corresponds to the bid-side. 
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Bid-side log(P) 




Figure 5: The time series plots of the deseasonalized bid- and ask-side liquidity 
factors, In fir and In /3 + , for TDC during January 2005 - March 2005. The lower 
picture gives the contour plots of the corresponding bivariate kernel-density 
estimate. The time interval is 10 minutes. 
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can be computed numerically by solving the initial value problem 



4 At = I + AAt 

using techniques of numerical integration^. The volatility matrix S can be 
obtained with a similar technique using the identity (0 and the estimated values 
of A and Vi 

We used the above procedure to fit Model to order book data with the 
sampling interval h equal to 10 minutes. Table [H gives the results for TDC 
and Moeller-Maersk. The estimates of the matrix A are consistent with the 
empirical autocovariance structures and the observed mean-reverting behavior 
of liquidity; see e.g. [5] or [16]. This shows up as the negative eigenvalues of A 
given in Panel B of Table [2j The negative eigenvalues imply, in particular, that 
the liquidity factors are mean reverting. In the words of Biais et al. [SJ page 
1657], "investors provide liquidity when it is valuable to the marketplace and 
consume liquidity when it is plentiful" . This can be explained e.g. by the so 
called "stimulated refill" effect which means that when liquidity is reduced by 
a market order on one side of the book new limit orders appear more frequent 
on that side; see e.g. [B] or [7]. 



4 In this study, we employ MATLAB routine ODE45. 
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The eigenvalue which is close to zero corresponds to the mid-price which 
is usually found to be nonstationary. In our three-month dataset, even the 
mid-price appears stationary in some cases. It can be seen as a consequence of 
the mean reversion property of the mid-price which is often reported for high- 
frequency data of stock prices; see e.g. [2] or Harris and Panchapagesan [T5] . 
This corresponds to the negative sign of the first element on the first row of A. 

The signs of the second and the third element on the first row of A suggest 
that the order book shape affects mid-price movements. In particular, the neg- 
ativity of the second element means that a steep bid-side lowers the drift of the 
mid-price; see Model[5J On the ask side, the effect is reversed. This corresponds 
to the so called "crowding out" effect which is supported by several studies of 
limit order markets; see e.g. Parlour [27], Griffiths et al. [18], Ranaldo [30] and 
Cont et al. [11] . It means that increased supply on one side of the book results 
in more aggressive incoming orders on that side and less aggressive orders on 
the opposite side. In our reduced form model, this shows up as increased mid- 
price drift when the ask-side liquidity, in terms of the liquidity factors, is high 
compared to the liquidity on the bid-side and vice versa. The results were found 
to be consistent across different stocks and different sampling intervals. 

4.4 Resiliency 

The observed mean-reversion properties of the order book models may be in- 
terpreted as a form of market resiliency. After a liquidity shock, the market 
mechanisms work by adjusting supply so that the liquidity properties of the 
book tend towards an equilibrium state. We will illustrate this with impulse 
response analysis where one of the variables of an estimated model is perturbed 
from its equilibrium value and then observing how the dynamic structure of the 
model operates to restore the equilibrium. 

As noted in the previous section, the estimated models are stationary which 
implies the existence of the limit £oo = lim^oo E£t. Taking expectations on 
both sides of (J2J we see that the expectation £t = E£t satisfies 

The equilibrium value £oo solves A£,+a = which has a unique solution when the 
eigenvalues of A are all zero. For TDC, we get ^ = (5.5069, -0.0084, -0.0229) 
and for Moeller-Maersk £«, = (10.9150,0.1148,-0.1735). Equilibrium proper- 
ties of limit order books have been studied also e.g. in Luckock [35] and Cont 
et al. Qj] but their results are not directly comparable to the above since the 
liquidity factors f3 + and /3~ have been deseasonalized; see Section |4"T2"1 

Figure [7] plots the impulse responses for TDC and Moeller-Maersk. In both 
cases we perturbed the liquidity factors /3 t ~ and /3 t + and tracked the effects on 
the mid-price drift and the liquidity factors. Positive deviations of the liquidity 
factors from their equilibrium values can be interpreted as liquidity shocks which 
reduce the liquidity on one side of the book. The first row of figures draws the 
responses to a shock on the bid-side and the second row on the ask-side. Each 
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Table 2: Parameter estimates for TDC and Moeller-Maersk. The panel A gives 
the parameter estimates for the linear equation <i£ = (At; + a)dt + T,dWt in 
Model[2j where £ = (In s, In /3~ , In j3 + ) . Panel B gives the spectral decomposition 
of the matrix A for both stocks. 

Panel A: Parameters for Model [2] 

Parameter TDC Mocllcr-Macrsk 



-0.0014 
-0.5132 
-0.9445 



-0.0003 
-0.2466 
-0.0133 



0.0004 
-0.0035 
-0.1952; 



0.0000 
-0.0001 
0.0000 



-0.0001 
0.1963 
0.0492 



0.0000^ 

0.0492 

0.1254/ 



'0.0080\ 
2.8240 
^5.1970/ 



-0.0002 
0.9247 
-1.3643 



-0.0001 
-0.6417 
0.1465 



0.0001 
0.1336 
-0.7363, 



0.0000 
0.0001 
-0.0002 



0.0001 
2.0846 
-0.0492 



-0.0002> 
-0.0492 
1.7570 



0.0022 
-9.9966 
14.7467, 



Panel B: Spectral decomposition of A 



TDC Moeller-Maersk 



-0.1893 0.0007 -0.0023\ /-0.4524 0.0002 0.0000^ 

Eigenvectors ( 0.3857 0.9664 -0.0424 -0.4984 0.5651 0.7995 

0.9030 0.2569 0.9991 / \ 0.7395 -0.8250 0.6006 y 

-0.0029\ /-0.0005^ 

Eigenvalues ( -0.2479 -0.8364 

-0.1925/ V -0-5413; 
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shock is defined as an increase of one standard deviation from the equilibrium 
value of the liquidity factor. The responses are given in terms of the development 
of the median and the 95%-confidence interval. 

A positive deviation in (3~ results in significant decreases in the subsequent 
drift of the mid-price. The decreased drift can be interpreted as decreased 
average returns. It takes about 30 minutes for /3 t ~ to return to its equilibrium 
value. There is practically no effect on the ask-side liquidity. The results are 
similar when perturbing ask-side liquidity factor j3^~ but the effect on the mid- 
price is reversed. This kind of behavior of the order book is consistent with the 
crowding out effect which, in microstructural analysis, is usually described in 
terms of incoming orders. 

5 Conclusion 

This paper proposes a reduced form approach for modeling limit order markets. 
The models are obtained by augmenting classical perfectly liquid market models 
by few additional risk factors that describe the liquidity properties of the order 
book. The proposed models are easy to calibrate and to analyze using standard 
techniques for multivariate stochastic processes. Despite the simple structure, 
the models are able to capture several empirically observed properties of order 
book dynamics. 

The modeling approach was illustrated by fitting three-dimensional continuous- 
time order book models to limit order data from Copenhagen stock exchange. 
The underlying risk factors in the models were the mid-price, and two liquid- 
ity factors describing the overall liquidity on each side of the order book. The 
estimation results support the chosen parameterization by revealing significant 
interactions between the liquidity factors and log-returns on the mid-price. The 
estimated model provides also a simple quantification of market resiliency with 
respect to liquidity shocks. 

While the proposed parsimonious approach will inevitably fail to capture 
some details of limit order markets, we believe that its simplicity will be an 
advantage when studying the costs of dynamic trading strategies. The mono- 
tonicity and convexity of the cost functions in the model are essential in that 
respect; see e.g. Almgren and Chriss [2J, Huberman and Stanzl [21] or Alfonsi, 
Schicd and Fruth [T]. Moreover, in situations where analytic solutions are not 
available, the low-dimensionality of our models makes them amenable to nu- 
merical techniques such as those proposed in Bertsimas and Lo [3] or Koivu 
and Pennanen 25J. This is currently under investigation and will be reported 
elsewhere. 
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(b) Moeller-Maersk 

Figure 7: Impulse responses for TDC and Moeller-Maersk order book models. 
The first row of figures gives the responses to an impulse on the bid-side impact 
factor. The response variables are the mid-price drift and the logarithmic bid- 
and ask-side liquidity factors. The second row gives the responses to an impulse 
on the ask-side liquidity factor. The time unit is 10 minutes 
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